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ABSTRACT 

We consider the description of cosmological dynamics from the onset of 
inflation by a perfect fluid whose parameters must be consistent with the 
strength of the enhanced quantum loop effects that can arise during infla- 
tion. The source of these effects must be non-local and a simple incarnation 
of it is studied both analytically and numerically. The resulting evolution 
stops inflation in a calculable amount of time and leads to an oscillatory uni- 
verse with a vanishing mean value for the curvature scalar and an oscillation 
frequency which we compute. 
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1 Introduction 



Although it is not yet known how to account for the late time acceleration of 
the universe [H [2] , it is by now quite clear that an adequate period of approx- 
imately exponential expansion - inflation [3] - provides a simple and natural 
explanation for the homogeneity and isotropy of the large-scale observable 
universe [3], and also is in satisfactory agreement with the primordial density 
perturbations spectrum [5J. This inflationary phase is usually realized by a 
scalar field but this is not a necessity. 

Consider the pure gravitational equations of motion: 

Gfj,y = Rf^u — -j^QfivR = —Ag^iv ■ (1) 

If A is assumed to be positive, the "no-hair" theorems imply that classi- 
cally the local geometry approaches the maximally symmetric solution at 
late times [B]; this solution is de Sitter spacetime and, thus, A-driven infla- 
tion is intrinsic to ([T]). 

As dl]) shows, classical gravitation without matter is a theory which only 
"knows" about the cosmological constant A ; Newton's constant G sets the 
strength of quantum effects. The corresponding mass scales are the Planck 
mass Mpi - associated with G - and the mass M - associated with A: 

We restrict ourselves to scales below the Planck mass so that the dimension- 
less coupling constant e = GA of the theory is small: d 

M < Mpi ^ e = GA < 1 . (3) 

The quantum behaviour of gravity in de Sitter {dS) spacetime has been 
studied in perturbation theory and in the infrared [ZIIH]- In particular, the ex- 
pansion rate H{t) decreases by an amount which becomes non-perturbatively 
large at late times: % 

deSitter =^ H{t) c:^ Ho [l - e^[#{Hot) + 0(1)] + 0{e^)} , (4) 



^Hellenic indices take on spacetime values while Latin indices take on space values. The 
Hubble constant is 3Hq = A. Our metric tensor has spacelike signature and our curvature 
tensor equals: i?"^^^ = r"^^_^ + r% T^^^ - (/z ^ i/). 

^This is a very mild restriction on the range of scales; for instance, if M ^ 10^^ GeV 
we get that e = GA - ^ - IQ-i^. 

^The Hubble parameter H{t) shall be defined in Section 2. 
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where 7^ is a positive pure number of 0(1). The underlying physical mech- 
anism is the production of infrared quanta out of the vacuum due to the 
accelerated expansion of spacetime. Such a production can only occur for 
particles that are light compared to the Hubble scale without classical confor- 
mal invariance; gravitons and massless minimally coupled scalars are unique 
in that respect. 

The factor of Hq t which appears in expression (4) is known as an infrared 
logarithm because it derives from infrared virtual particles and because Ho t 
is the logarithm of the de Sitter scale factor. Any quantum field theory which 
involves undifferentiated gravitons or massless minimally coupled scalars will 
show infrared logarithms in some correlators at some order in the loop ex- 
pansion. If the interaction contains undifferentiated gravitons or massless 
minimally coupled scalars, along with any number of other fields, then each 
new factor of the couphng constant squared can produce at most additional 
infrared logarithms. For example, the fundamental interaction of quantum 
gravity in de Sitter background takes the generic form [9]: 



where h^j^^, is the fluctuating graviton field. Thus, one can get at most one 
extra infrared logarithm for each additional power of G [S]. 

The operator under study also has an effect. For example, because there 
are two derivatives in the invariant measure of acceleration [10] whose expec- 
tation value gave expression (4), the general form of such corrections is: 



where i stands for the loop order and where the constants cik are pure num- 
bers of 0(1). 

Because e is constant, whereas HqI grows without bound, infrared loga- 
rithms eventually lead to a breakdown of perturbation theory. In quantum 
gravity on de Sitter background this occurs after about e-foldings. To 
evolve further requires a non-perturbative technique such as summing the 
series of leading infrared logarithms: 




(5) 
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Starobinskii [TT] has developed a stochastic technique which exactly repro- 
duces the leading infrared logarithms of scalar potential models [H] , and which 
can be used to sum them whenever the scalar potential is bounded below 
^12j . Starobinskii's technique has recently been extended to include models 
in which the scalar interacts with other fields such as a Yukawa fermion [13] or 
electromagentism [T3] . The late time limits of the vacuum energies of these 
scalar models exhibit a broad range of possibilities for what the quantum 
gravitational sum might give: 

• Scalar potential models which are bounded below show a small, con- 
stant increase of the vacuum energy IT5] : 

• Scalar quantum electrodynamics experiences a small, constant decrease 
of the vacuum energy [15]; and 

• Yukawa theory engenders a decrease of the vacuum energy which grows 
without bound [13]. 

We would ultimately like to compute how quantum gravity affects late 
time cosmology by employing Starobinskii's technique to sum the series of 
leading infrared logarithms. There has been some progress in this area [16] 
but the full solution is not yet in sight. A more modest approach is to antic- 
ipate the solution by attempting to guess the most cosmologically significant 
part of the effective field equations guided by our understanding of the per- 
turbative regime at leading logarithm order. The resulting model could be 
regarded as a worthy object of study in its own right, just as one views the 
many classes of scalar-driven inflation models, without feeling any need to 
derive them from fundamental theory. It may even be that our study will un- 
cover some general feature of any successful model which can, it turn, guide 
the fundamental derivation. 

It should be noted that doubts have been raised about the possibility 
of any infrared contribution to the quantum gravitational vacuum energy 
[T71 [T8] . However, these doublts are difficult to reconcile with the fact 
that scalar models certainly show infrared corrections to the vacuum energy 
[T^ IT^ [m |T3], and with Weinberg's observation that infrared logarithms 
contaminate the power spectrum of scalar-driven inflation [IQj. General the- 
oretical arguments have also been advanced to show that de Sitter must be 
unstable in quantum gravity [20] . 

One feature which complicates evaluation of these arguments is the in- 
tractibility of quantum gravity at any order, and the fact that the onset 
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of this particular effect occurs at two loops. The latter fact must be so 
because screening represents the gravitational attraction between virtual in- 
frared gravitons which have been ripped from the vacuum. The production 
process is a one- loop effect, so the gravitational response to it cannot occur 
until the next loop order. Three separate computations of the graviton 1- 
point function have confirmed that there is no one- loop effect [2T], [221 123] , 
and the same conclusion can be reached from taking the de Sitter limit of 
scalar-driven inflation [25] . This difficulty of performing explicit com- 
putations is one more reason why it might be desirable to study quantum 
gravitational screening from the perspective of the effective field equations. 

In the present paper we shall use the physical principles responsible for 
the non-trivial quantum gravitational back-reaction on inflation to construct 
a phenomenological model which we can then directly evolve. Therefore, 
we wish to construct an appropriate effective conserved stress-energy tensor 
Tij.u[g] which will modify the gravitational equations of motion in the 
usual way: 

C^. = -A^,.. + S7TGT^,[g] . (8) 

Our stress-energy tensor must be a non-local functional of the metric as 
dictated by the nature of the effect we wish to describe. It can be conveniently 
parametrized as a "perfect fluid" : 

T^^^a] = {.P + p)u^,Uy + pg^y , (9) 

so that to completely determine it we need the following three ingredients: 

(i) the energy density p as a functional of the metric tensor p[g\{x), 

(a) the pressure jo as a functional of the metric tensor p[g]{x), 

(in) the 4- velocity field 11^ clS Qb functional of the metric tensor u^[g]{x), 

chosen to be timelike and normalized: 

g^'''u^u, = -l =^ u^'u^,, = Q . (10) 

Section 2 describes a simple ansatz for the effective stress-energy tensor 
T^v[g] and the set of equations it leads to for homogeneous and isotropic 

''Note that the slow roU suppression one finds for corrections to the background in this 
Hmit merely means that some of the fields must be differentiated, as in the VG h dh dh 
vertex. Hence, one can only get a single extra infrared logarithm for each extra G, as 
opposed to the three infrared logarithms that would be possible if the vertex had been 
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spacetimes. Section 3 presents numerical results obtained by discretizing the 
relevant evolution equation. To the degree that is possible, the dynamics 
of homogeneous and isotropic evolution in the presence of T^y[g\ is studied 
analytically in Section 4 and - wherever a comparison can be made - the 
excellent agreement with our numerical study is noticed. Section 5 discusses 
late time evolution and the possible modifications it implies on our ansatz. 
Our conclusions comprise Section 6. 

2 A Physical Ansatz 

It is clearly impossible to uniquely fix the functional form of the effective 
stress-energy tensor solely from the physical requirements and correspon- 
dence limits we have at our disposal; only quantum field theory could, in 
principle, provide such an answer. However, we can try to obtain a simple 
ansatz and then analyze its implications. 

• Why a Perfect Fluid? 

One might think that a good way of discussing quantum corrections to the 
field equations would be in terms of an ansatz for the effective action. How- 
ever, the "in-out" effective action - being non-local - does not give causal 
effective field equations we can use to study cosmological evolution. The "in- 
in" effective action of the Schwinger-Keldysh formalism does produce causal 
effective field equations, but this results from subtle cancellations between 
different off-shell fields, making it difficult to identify promising candidates 
for the off-shell effective action. For a class of simple non-local effective ac- 
tions, it is possible - using a partial integration "trick" - to obtain causal 
effective field equations [27j but it is impossible to restrict the non-local effects 
to the past. The chain of action and reaction that is part of any single-field 
variational formalism implies that screening the cosmological constant by an 
inverse differential operator acting on some curvature scalar will inevitably 
lead to a variation of that curvature scalar appearing at the current time 
in the effective field equations. Hence one always gets a renormalization of 
the effective Newton constant. We wish to avoid this, so we specify the ef- 
fective field equations directly and insist that the non-local screening of the 
cosmological constant both remains in the distant past and also does not 
change the Einstein tensor so that Newton's constant is not renormalized. 
The problem with this procedure is that we must enforce conservation. We 
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selected the perfect fluid form for the effective stress tensor both because en- 
forcing conservation is straightforward and because our perturbative studies 
of quantum gravitational screening indicate that this form suffices to capture 
the leading infrared logarithms we seek to reproduce [26] . 

• Implications of Conservation 

The "perfect fluid" parametrization Q of T^u[g] allows us to completely 
determine it from the three quantities it contains: the scalars p, p and the 
4- vector m^. Because of the normalization (fTOj) . only three of the components 
of are algebraically independent. Thus, Tf^u[g] contains five independent 
quantities in total. Conservation provides four equations and allows us to 
determine any four of these quantities in terms of any one. It turns out to 
be more convenient to specify the induced pressure functional p[g] and then 
use conservation to obtain the form of the induced energy density p[g] and 
4- velocity Ufj_[g] up to their initial value data. 
The fundamental equation: 



D^T^. = , 



implies: 

dyp + Uy{u - D + D ■ u){p + p) + {p + p){u- Duy 
By contracting into (fT2!) and using (fTOl) . we get: 



u ■ dp 



{p + p)u^ 



u-dp = —{D-u){p + p) 



'111 



(12) 



(13) 
(14) 



Then, by substituting (HM into the conservation equation flT2|) . the following 
equation emerges: 



Equations 



{p + p) u- Duy = -{dy + UyU- d)p 
can be used to accomplish our goal. 



(15) 



• Requirements on the Pressure 

(i) The initial value requirement 
Our gravitationally induced source should not disturb the basic nature of 
the pure gravitational equations ([T]). The latter can be evolved from the 
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initial spacelike surface knowing only the metric and its first time derivative. 
This property of gravity must be retained in the presence of the source and 
constrains both the local and non-local parts of its functional form |28] ; for 
instance, any local parts in T^jj[g\ can contain at most second time deriva- 
tives of the metric. 

(a) The non-locality requirement 
We argued that the physical effect responsible for gravitationally inducing 
Tp.v[g\ is inherently non-local and, therefore, our source must be non-local. It 
is important to mention that this conclusion can also be reached by noting 
that no local modification of pure gravity can prevent de Sitter spacetime 
from being a solution of the field equations eternally [2Sj. Any local modifi- 
cation simply changes the initial Hubble constant ifo and can be absorbed by 
the cosmological constant counterterm 5K to leave no change and de Sitter 
spacetime as a solution for all time. Thus, the important part of the induced 
stress-energy tensor must be non-local. 

(Hi) The simplicity requirement 
A simple non-local operator at our disposal is the inverse of the scalar 



and a simple scalar it can act on is the curvature scalar R. Hence, we shall 
explore ansatze in which the pressure is a function of the quantity X[g\: 



(iv) The correspondence requirement 
Our gravitationally induced source should reproduce the perturbative results 
obtained in de Sitter spacetime. 

• Cosmological Spacetimes 

The large-scale homogeneity and isotropy of the universe selects Friedman- 
Robertson- Walker {FRW) spacetimes as those of primary cosmological in- 
terest; their line element for zero spatial curvature equals in co-moving coor- 
dinates: 



d'Alembertian: 




(16) 



ds^ = g^u{t) dx^dx" = —dt^ + a^{t) dx ■ dx . 



(18) 



^Our scalar d'Alembertian is defined with retarded boundary conditions. 
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Derivatives of the scale factor a{t) give the Hubble parameter H{t) - a mea- 
sure of the cosmic expansion rate - and the deceleration parameter q{t) - a 
measure of the cosmic acceleration: 



Hit) 



a{t) 



ait) 
a{t) d{t) 



d 



-1 



Hit) 



(19) 
(20) 



d2(t) H^{t) ' 

For these spacetimes the stress-energy tensor iQ takes the form: 

uouoip + p) - p = p , (21) 
, (22) 



^00 

T- 



UiUj (p + p) + QijP = gijP 



(23) 



An immediate consequence of isotropy and the normalization condition ([TUj) 
is: 



5.' 



The Ricci tensor and Ricci scalar become, respectively: 

3a 



and: 



Roi 
Rij 

R 



a 



3H' 



3H 



, 
d 
a 



+ 



2d 



;,2 



9ij 



6d ^ 6(2^ 



a 



(l2//2 



QH 



(24) 

(25) 
(26) 

(27) 
(28) 



8]), the non-trivial gravitational equations of motion 

(29) 
(30) 

(31) 



In view of 
([HD take the form: 

3//2 = A + SnGp , 
-2H-3H^ = -A + 8nGp 

while the conservation equation ffTTl) becomes: |^ 

p = -3Hip + p) . 



^The FRW conservation equation ([3T|) can also be derived directly from the equations 
of motion (l29l[30l). 
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The latter implies that: 



p{t) = -p{t) + 



dt' a\t') Pit') . 



(32) 



a\t) Jo 

When acting on functions which only depend on co-moving time, the 
scalar d'Alembertian (fT6l) for FRW geometries equals: 

□ = -[d^ + 3Hdt) , (33) 

so that its inverse is: 

- = - / rft' 



dt" a^t'" 



lo a^{t') Jo 
Consequently, the source X can be written as follows: 



(34) 



X 



R 



* 1 
dt' 



a3(t') Jo 



dt" a%t") l2H\t") + QH'\t") . (35) 



Note that we have taken the initial time to be at t = 0. 



• The de Sitter Correspondence Limit 

If we define inflation as positive expansion, i.e. H{t) > 0, with negative 
decelaration, i.e. q(t) < 0, a locally de Sitter geometry provides the simplest 
paradigm. It is characterized by constant Hubble and deceleration parame- 
ters, and a scale factor of a simple exponential form: 



Qds{t) 



a-dsit) 



(36) 



The source X^s can be computed by first obtaining the inverse d'Alembertian 
and curvature scalar from and (1251) respectively: 



1 /■* , /■*' 

- = - dt' e-3^°* / dt" 

□ dS Jo Jo 

Rds = , 
and then acting the former on the latter: 



,3Hot" 



(37) 
(38) 



X. 



dS 



□ 



dS 



Jo Jo 



-4//o 
-AHot 



t dt' 




lo 





1 -e" 



-3Hot 



(39) 
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For large observation times we get: 



-41n[a,5(t)] + 0(1) , \n[ads{t)] = Hot . (40) 



/^From fl40l) we already see that our source X can give the infrared logarithm 
in our perturbative result (jlj). We should mention at this stage that explicit 
perturbative computations for a massless minimally coupled self-interacting 
scalar field give the following leading logarithm result ^29] : 



e-i 



leading log =^ H\t) = - \ l-GAj2 GA \n[ads{t) ] ^ ^ \ (41) 

3 I e=2 



where i is the loop order and are pure numbers. Now consider the equation 
of motion ( !30l) for de Sitter spacetime: 

- 3H^ ~ + 8nGp[gas] ■ (42) 

The leading logarithm form of the induced pressure immediately follows: 

leading log =^ p{t) = ^ hi GK\n[ads{t)] ■ (43) 

e=2 

In view of the above, the following ansatz for the gravitationally induced 
pressure p[gds\if) reproduces (jl]) up to a numerical coefficient: 

p[gdsm = A'f[-eXds]{t) = A' f[-{GA){n-'RU]{t) , (44) 

where / is some monotonically unbounded function. The reason for requiring 
/ to be unbounded is dictated by the similar behaviour seen in our lowest 
order perturbative result (jlj). Moreover, since the source X^s is increasingly 
negative definite as can be seen from fl39|) . any function / satisfying: 



f[-eXds] = -eXds + 0{e') , (45) 
will be increasingly positive definite. 



^According to our perturbative result ^ , Hds (t) is subdominant since its time deriva- 
tive must eliminate one infrared logarithm ln[a(i5(t)] without affecting the corresponding 
factor of e = GA. 
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Finally, successively apply equations (HSl) and (HUj) to the equation 
of motion 



~ H^{l - Svre [ -e + 0(e2) ] } 

^Hf{l - 32ne^ \n{ads) + ^(e^) } . (46) 

Up to a positive numerical coefficient, |fl the requisite agreement is achieved. 

• The General Ansatz 

Our physical requirements and correspondence limits have led us to the fol- 
lowing ansatz for the gravitationally induced pressure p[6'](2;) in a general 
geometry: 

p[g]{x) = A'f[-eX]{x) , ^ = 5^ > (47) 
where the function / satisfies: 

f[-eX] = -eX + 0{e') . (48) 

To completely determine the induced stress-energy tensor T^^[g]{x) we need 
the energy density p[g]{x) and the 4- velocity u^[g]{x). Given the pressure 
]9[g'](x), we can obtain the other two quantities via stress-energy conservation 
up to their initial value data. 

An explicit cosmological model needs an explicit function /. Out of the 
plethora of functions satisfying fBSl) we shall select simple ones and define: 

(z) the linear model : f[—eX] = — eA , (49) 
(ii) the exponential model : /[— eA] = e^'^^ — 1 . (50) 

It is the dynamical evolution of the explicit cosmological models that will 
determine whether, as we expect, they will naturally stop inflation. After- 
wards, the induced source should "turn-off" and the universe should enter a 
radiation dominated epoch. 



®In the pure number /i2 was calculated and was found to equal j^^- Since this 
was done within the context of a simplified quantum gravitational theory, we shall instead 
assume it equals one and, therefore, it does not alter the coefhcient of the —eXds term in 
(HI). 
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3 Numerical Results 



Because the non-local, non-linear equations we have proposed are too com- 
plicated to solve exactly, we shall evolve them numerically. For that purpose, 
it is preferable to use the {ij) equation of motion due to its linearity in the 
highest derivative: 

1 



2H + 3H^ = 3H^{l - 87V€f[-eX]} , X = -^R , (51) 

where, as we have already mentioned, Hq is the Hubble parameter at the 
onset of inflation and e = GA = 3GHq is the dimensionless coupling constant 
of the theory. 

The discretization of (15T|) involves: 

(i) Constants, 

step size in Hubble units =^ S = HoAt 
coupling constant =^ e = GA , 

(ii) The basic variables and initial value data,^ 

a{t) — > a{i At) = e^' , &o = , 



(52) 
(53) 



(54) 



{p + p){t) 



[p + p]{iAt) 



[p + P 







[P + P]^^, = e-''^[[p + pl-eAX,f[-eX,]] . (55) 

All quantities of interest as well as their initial values can be determined from 
dSlES]): 

(i) Dynamical quantities, 

Ah 



m — 


H{iAt) 


m — 


H{iAt) 


m — 


* q{iAt) 


Rit) — 


^ R{i At) 


X{t) — 


X{iAt) 







At 
A\ 
At 2 



Ah = 
A%, 



1 
1 



A'h 



At2 



Abf ' 

QA\ + 12(A6i)2 
Xi + AX, , 



h+i - h , (56) 
E Ah+i - Ah , (57) 
(58) 
(59) 



AX, - 12 (Ah)^ - 6 A%i 



(60) 



^It is of course equivalent to consider (p + p) instead of p because of stress-energy 
conservation; see equation ()31|) . 
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(a) Initial value data, 



Abo 



= 6 
AXo 



A\ = 







The discretized evolution equation fl^ reads: 



A'%, 



3r 
2 



(61) 
(62) 



(63) 



It was numerically integrated using Mathematica for the following choice of 
the input parameters and of the step range: 



1 



1000 



1 
200 



i e [0, 350000] 



(64) 



Moreover, for the function / we chose the one corresponding to the exponen- 
tial model fl50D: 



e"- 1 



f-\x) = ln(l + x) 



(65) 



The associated critical point Xc 
uj - defined in (175]) - are: 



-eXcr - defined in (!70|) - and frequency 



UJ 



e5_ 
At 



eS_ 
At 



- In (l + ■ 

72n (l + 



1 



57re 
1 

8^ 



438.50 , 

2.25 X 10"^ 
At 



(66) 
(67) 



Our results are presented in the set of graphs that can be found in the 
very end. Some comments are in order: 

• After the onset of and during the era of inflation, the source X{t) grows 
while the curvature scalar R{t) and Hubble parameter H{t) decrease. 

• Inflation ends and the time when this occurs is the time when the decel- 
eration parameter q{t) goes from negative to positive values. 

• During the era of oscillations: 

(i) The oscillations of R{t) are centered around R = 0, have an envelope 
behaving like t~^ and a frequency u in agreement with fl67|) 
(a) Although there is net expansion, the oscillations of H{t) take it to small 
negative values for small time intervals. The presence of these short deflation 
periods is a novel feature of the model and may have consequences on the 
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primordial perturbation spectrum. 

(Hi) The oscillations of H{t) show that there is almost no difference between 
R{t) and 6H and, therefore, the term in R{t) proportional to H'^{t) is in- 
significant during this era. 

(iv) The oscillations of the scale factor a{t) are centered around a linear 
increase with time. 

It is important to note at this stage the excellent agreement of the ana- 
lytical results derived in Section 4 with their numerical equivalents. The two 
basic parameters to concentrate are, 

• Criticality: It occurs at step i = 160942, as the detailed data of Figure 12 
indicates, and at that point: 

X[160942] = -438.50 , ^[160942] = 0.50 , (68) 

in complete agreement with the analytical predictions fl66p and 0911) respec- 
tively. 

• Oscillation frequency: /^From the detailed data of Figure 4 we conclude 
that six oscillations have occured between steps i = 174291 and i = 342478. 
Hence, we have: 



342478 - 174291 , 2n 

T = At = — 

6 io 



2.24 X 10"^ , , 
"=^^> (69) 



which compares very well with the analytical prediction ( 1671) . 



4 Analytical Results 

With the evolution equation ( 15T1) as a starting point, we can analytically 
derive some results for the physical system under study. Our ansatz restricts 
the function / to be monotonically unbounded. Therefore, there exists a 
critical point Xcr such that: 

1 - 8nef[-eX,r] = =^ X,, = . (70) 

Inflationary evolution dominates roughly until we reach the critical point. 
Close to the critical point the induced pressure p is small and, thus, it makes 
sense to expand / around its critical point and use the resulting perturbation 
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theory for the subsequent evolution: 

2H + 3i/2 = 2,Hl { 1 - Svre /[-eX,, - e(X - X^,)] } 

= 3//o'|l - 87re(/[-eX,,] - e{X - X„) f[-e X,, 

+ 0(e\X-X,^,f) I 

- 2A'Ke^Hl{X-X,r)n-eX„] . (71) 

Neglecting all higher order terms is a superb approximation given the very 
small realistic values of e. 

Moreover, using (128|1 we rewrite the co-moving time derivative of the 
Hubble parameter as: 

H = - 2H'^ . (72) 



Consequently, the evolution equation flTTl) becomes: 

-R + 3H^ ^ -72n {tH^f {X - X,,) , (73) 
where we have defined: 

. . -\^n-,X]\^__^^^ . (74) 

Action of the d'Alembertian operator fl33l) on (173!) gives: 

R + 2H R + {u"^ - H)R + [3H^R - 3QH^] ~ , (75) 
with the understanding that: 



u = e//oV72vr/^, • (76) 

It will turn out that the term in brackets is subdominant and we can focus 
our attention to the differential equation: 

R + 2HR + - i/) - , (77) 

which describes a damped oscillator. To solve the above equation we first 
scale out the Hubble friction term by defining: 

R = -S , (78) 
a 
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so that (1771) becomes: 

S + {u^ -2H - H^) S c:^ , (79) 
and, in the large time limit, is solved by: 

S{t) ~ Kism{iut + ip) , :$> \ - 2H - . (80) 

Hence, by using (!78|) and fl28|) . we have: 



sin(cjt + (/?) aR = Q^{Ha) + 6H^a ~ Q-^{Ha) . (81) 



/^From (18111 we immediately conclude: 



K 

a{t) K2 - cos{ut + ip) , (82) 
o u 



a{t) ~ + fsTs t - ^ sin(a;t + ip) . (83) 



6 uj 

and our large time results are: 



N 1 Ki cos( u>t + Lo) ^/1\ . 

lim H(t) ^ — + O(-) , 84 

t»i ^ ^ t QK2UJ t ' ^ ' 

xV/ N Ki sin(u;t + 09) ^/1\ /^^s 

^1 !!!lM±^ + o(i) ^ (86) 

i>l ^ ' K2 t H^^ ^ ' 

iFiom. the asymptotic solution that we just obtained, the physical picture 
that emerges so far is that of a universe in which, as the exit from the infla- 
tionary era approaches, oscillations in R become significant] their frequency 
to is given by fl76l) and, according to (1861) . their envelope is linearly falling 
with time. 

We can get further insight by computing the deceleration parameter q at 
criticality. This is most easily done by considering the equations of motion 
f l29ll30|) from which we deduce that: 

-2H = 8ttG{p + p) = A87re^^ . (87) 



^"These results justify our ignoring the bracketed term in equation ([75|l and the H^a 
term in equation ()8ip . 
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Then, we rewrite f l29p in a form convenient for our purpose: 



3i/^ = A + ^TiGp 



A 



A 



1 - 8^^A^ + S^^^ 



(88) 



and use fl87H88l) to express the decelaration parameter (l20l) in the following 
way: 



H , 3 



A87re£^ 



8ne- 



+ Sire 



p+p 



A2 ^ A2 

The definitions (1701) of criticality and fHTl) of pressure imply: 

1 - 87re/[-eX,,] = 1 
and allow us to conclude that: 



A2 



-I 



1 

+ 2 



(89) 



(90) 



(91) 



At the onset of inflation go = Since by the time the universe arrived 
at the critical point the decelaration parameter had already reached positive 
values - Qcr = +^ - the epoch of inflation ended before the universe evolved 
to the critical time. 

To investigate whether the dynamical system is underdamped at the crit- 
ical point, we isolate all terms of the full equation (1751) that affect the fre- 
quency: 



frequency terms 
We have already seen that: 



(^uj'^ - H + 3H^) R 



■1 



Hr 



+ 



1 



HI ■ 2 

leading to a positive frequency determining coefficient: 

9 



+ - Hi > 



(92) 



(93) 



(94) 



Hence, at criticality the system is underdamped implying again that oscilla- 
tions start around the end of inflation. 
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It is also interesting to work out an approximate but direct relation be- 
tween the frequency uj and the Hubble parameter H at the critical time. The 
starting point is equation fl87|) and the approximation consists of evaluating 
its right hand side in de Sitter spacetime: 

-m ^ ^(P + ^^)L = 32vr (ei/o) , (95) 
and afterwards at criticality: 

H,, ~ -167r(e//o)Vd5 • (96) 
Comparison of (!96|) with (!76|) gives us the desired approximate expressions: 

a;2 27 ' 27 ' ^ ' 

where we have used (1931) to arrive at the second relation. Even with the 
eternal de Sitter assumption - which ignores the dimunition of the Hubble 
parameter H with time - the frequency uj is the larger quantity; in the real- 
istic case the ratios fl97|) would be much smaller. What we can conclude is 
that, in our evolution, the inequality < \H\ < u"^ is well justified. 

-k Identities 

In addition to the equations presented throughout the main text, various 
of the following expressions have been used to obtain the results of this 
subsection: 



R 



^ S-HS 



R = -\S-2HS + H^S-HS 



(98) 



-4Hot , X^s - -4i/o , (99) 
1 . _ 1 



5 After Inflation 

The homogeneous and isotropic evolution described in the previous two Sec- 
tions does not give a completely satisfactory end to inflation. The oscillations 
are no problem, but the average expansion a{t) ~ t is unacceptably rapid. 
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At that rate there would be no reheating and the late time universe would 
be cold and empty. Nonetheless, the same is true for scalar-driven inflation if 
one ignores the possibility for energy to flow from the inflaton into ordinary 
matter. We believe that energy will flow from the gravitational sector of our 
model into ordinary matter to create a radiation-dominated universe, just 
as it is thought to do for scalar-driven inflation. In that case, one should 
think of the total stress-energy as consisting of our quantum gravitational 
perfect fluid plus the energy density and pressure of radiation, with the latter 
described just as in conventional cosmology. 

An amazing possibility arises if this process can be shown to occur: our 
quantum gravitational correction cancels the bare cosmological constant and 
then becomes dormant during the epoch of radiation domination. To see this, 
suppose the deceleration parameter has the pure radiation value of q{t) = +1 
for times t > tr- In that case, the Hubble parameter and scale factor are: 

^ ^ 1 + 1H,f{t — t,f^ ^ ^ 

1 

a{t) = ar[l + 2Hoit - tr)Y , (102) 

where Hr and are their values at t = tr- During this phase the Ricci scalar 
is zero: 

g = +1 =^ R = QH+12H^ = . (103) 

Our simple source X{t) obeys the differential equation HX = R, so for t > tr 
it must be a linear combination of its two homogeneous solutions: 



Vt > tr =^ UX = 



X{t) = Xr + XrJ^dt' 



a 



r 



X 



Xr-jf I . (104) 

Hr ■ll + 2Hr{t-tr) 



ait') 

The only solution consistent with g = +1 is: 

Xr = X^r , Xr = ■ (105) 

The system is predisposed to reach nearly X^r in any case but one might 
doubt that evolution would enforce the exact vanishing of the second solution 



^^The equation CX ~ R remains true even if the total stress-energy includes radiation 
and/or matter contributions. 
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implied by = 0. However, note that the vanishing of the J^dt'a~^{t') 
homogeneous solution is attained by the purely gravitational evolution of 
Sections 3 and 4, which does not include energy transfer to matter. In that 
case i? 7^ 0, so there is a homogeneous contribution as well, but the fact that 
the oscillations are about Xcr means that the second homogeneous solution 
is completely absent. 

Having X{t) approach X^r within the context of a hot, radiation domi- 
nated universe would be a great success for our model, but the eventual tran- 
sition to matter domination poses problems. The onset of matter domination 
is really a gradual process but let us simplify the exposition by considering 
a sudden change from g=+ltog = +|at some time ty. During this 

matter dominated epoch the Hubble parameter and scale factor are: 



1 

+ 2 



H{t) 
a{t) 



Hr, 



1 + 2^m{t — tr. 



1 + -^Hjnit — t„ 



(106) 
(107) 



where Hm and H{tm) and a{tm), respectively, computed from the 

radiation dominated geometry fll01H102|) . During matter domination the 
Ricci scalar is nonzero: 



1 



+ 



R = 6H + 12H' 



(108) 



The resulting change in the source X{t) is: 
q - + ^ =^ 



AX{t) = X{t) - Xc 



3 



1 + -Hmit — tn 



0(1) . (109) 



To understand what is wrong with the change fll09p caused by matter 
domination, it is useful to recall our ansatz ( l47l) for the quantum gravita- 
tionally induced pressure: 

p[g]{x) = A^f[-GXX]{x) . (110) 

In the context of this ansatz there are two major problems with fll09p : 
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• The sign problem. It derives from the function f{x) in flllOp be- 
ing monotonically increasing and unbounded. Hence, pushing X{t) below 
Xcr -C results in positive total pressure, whereas observation implies nega- 
tive pressure during the current epoch [H |2]. Note that we cannot alter this 
feature of f{x) without sacrificing the very desirable ability of the model to 
cancel an arbitrary bare cosmological constant. 

• The magnitude problem. In one sentence, the magnitude of the total 
pressure produced by (11091) is vastly too large. The problem arises from the 
factors of the bare cosmological constant A in our ansatz (IllOp . The total 
pressure ptot is the sum of the classical contribution and our ansatz (IllOl) : 

Ptot = - ^|l - 87rG'A/[-G'A(X„, + AX)]| (111) 

^ - ^x(GA)V;AX . (112) 

Note that we need not include in ptot an additional contribution because non- 
relativistic matter has zero pressure. Comparing with the currently observed 
value jOnow of the pressure: 

Pnow — ~ o ^ -^now ; (H^) 



gives: 

Ptot fGAHo 



rr ! fcr^X ^ 10^^ X X AX , (114) 

where we have assumed Hq ~ 10^'^ GeV and i/now ~ 10"^^ eV. The deriva- 
tive /^^ is unity for the linear model (jinj) and of order {Gh)~^ ~ 10^^ for 
the exponential model fISUl) . so we expect /^^ to be at least of order one and 
possibly much greater. 

There is no way of addressing either problem without generalizing our 



ansatz (11 101) for the pressure. This necessarily takes us away from what can 
be motivated by explicit computation during the de Sitter regime. Although 
these issues will be analyzed elsewhere pQ], we shall mention the basic prin- 
ciples: 

• The magnitude problem arises because the constant A in (IllOp is about 
the square of the inflationary Hubble parameter rather than its late time de- 
scendant that could be 55 orders of magnitude smaller. Solving the problem 
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entails replacing one of these factors of A present in flllOp by some dynamical 
scalar quantity that changes as time evolves in a way that also preserves the 
original relaxation mechanism. 

• The sign problem arises because the Ricci scalar is positive during both 
inflation and matter domination. Again solving the problem involves a dy- 
namical scalar quantity that changes sign from inflation to matter domination 
and is still zero during radiation domination. 

6 Epilogue 

We have presented a simple ansatz for the most cosmologically significant 
part of the effective field equations of quantum gravity with a positive cos- 
mological constant. The quantum correction to these equations consists of 
a "perfect fluid" stress-energy tensor in which the pressure is specified as 
a non-local functional of the metric, and the associated energy density and 
timelike 4- velocity are determined by conservation. On the basis of simplicity, 
as well as correspondence with perturbative results in de Sitter background, 
we proposed that the pressure takes the form: 

p[g]ix) = A^f[-GAn-'R]ix) , (115) 

where □ is the scalar d'Alembertian and its inverse is defined with retarded 
boundary conditions. 

We studied homogeneous and isotropic evolution in this model, both nu- 
merically (Section 3) and analytically (Section 4). As long as the function 
f{x) is monotonically increasing and unbounded the qualitative behavior is 
the same: 

• Inflation is nearly de Sitter for a calculable period; and then 

• The Ricci scalar oscillates about zero with a calculable constant period 
and an amplitude that falls off like t~^. 

The universality of this behaviour was checked numerically by evolving func- 
tions f{x) all the way from linear to exponential, and Section 4 presents a 
derivation from the effective field equations. 

Of course (single) scalar-driven inflation contains a free function, the 
scalar potential V{ip), which can be fine-tuned to support a wide variety 
of expansion histories a{t). However, the generic evolution of our model in- 
volves two distinct features which single-scalar inflation can never reproduce: 

• During the oscillatory phase, the Hubble parameter H{t) actually drops 
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below zero for brief periods; and 

• The derivative of the Hubble parameter H{t) is positive for about half of 
the time during the phase of oscillations. 

The first feature is conducive to rapid reheating, while the violation of the 
weak energy condition implicit in the second is the hallmark of a quantum 
effect 

Prominent among the list of topics for future work is perturbations. We 
need to show that the dynamical scalar mode of our model releases the energy 
of oscillations into matter to reheat the universe. If this happens, then the 
quantum gravity sector will go quiescent during a long epoch of conventional 
radiation domination. The subsequent transition to matter domination might 
even give rise to something like the current phase of acceleration, but this 
requires modifications of our simple ansatz which shall be described elsewhere 

m- 

Another important topic for future work is to derive and solve the equa- 
tion for scalar perturbations, at least enough to compute the scalar power 
spectrum. One also needs that there be no long-range scalar force at late 
times. Moreover, the equation for tensor perturbations is unchanged by our 
"perfect fluid" model. We need only use the expansion history a{t) predicted 
by our model in order to compute the tensor power spectrum. 
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Figure 1: The evolution of the source X{t) over the full range for the exponential model. 
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Figure 2: The evolution of the source X{t) during the oscillatory regime for the 
exponential model. 
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Figure 3: The evolution of the curvature scalar R{t) over the full range for the 
exponential model. 
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Figure 4: The evolution of the curvature scalar R(t) during the oscillatory regime 
for the exponential model. 
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Figure 5: The evolution of the Hubble parameter H{t) over the full range for the 
exponential model. 
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Fi gure 7; The evolution of H(t) over the fuU range for the exponential model. 




Figure 8: The evolution of H{t) during the oscillatory regime for the exponential model. 
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Figure 9: The evolution of the deceleration parameter q{t) during the oscillatory regime 
for the exponential model. 




Figure 10: The evolution of the decelaration parameter q{t) around the end of inflation 
for the exponential model. 
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Figure 11: The evolution of the scale factor ratio [a(t)/a(150000)] during the oscillatory 
regime for the exponential model versus a linear interpolation. 
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Figure 12: Determining the critical point 1 — 87re/[— eX^] = for the exponential model. 



